Learning high-dimensional directed acyclic graphs with latent and selection variables

We consider the problem of learning causal information between random variables in directed acyclic graphs (DAGs) when allowing arbitrarily many latent and selection variables. The FCI (Fast Causal Inference) algorithm has been explicitly designed to infer conditional independence and causal information in such settings. However, FCI is computationally infeasible for large graphs. We therefore propose the new RFCI algorithm, which is much faster than FCI. In some situations the output of RFCI is slightly less informative, in particular with respect to conditional independence information. However, we prove that any causal information in the output of RFCI is correct in the asymptotic limit. We also define a class of graphs on which the outputs of FCI and RFCI are identical. We prove consistency of FCI and RFCI in sparse high-dimensional settings, and demonstrate in simulations that the estimation performances of the algorithms are very similar. All software is implemented in the R-package pcalg.Comment: Published in at http://dx.doi.org/10.1214/11-AOS940 the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org

Community based trial of home blood pressure monitoring with nurse-led telephone support in patients with stroke or transient ischaemic attack recently discharged from hospital.

BACKGROUND: High blood pressure in patients with stroke increases the risk of recurrence but management in the community is often inadequate. Home blood pressure monitoring may increase patients' involvement in their care, increase compliance, and reduce the need for patients to attend their General Practitioner if blood pressure is adequately controlled. However the value of home monitoring to improve blood pressure control is unclear. In particular its use has not been evaluated in stroke patients in whom neurological and cognitive ability may present unique challenges. DESIGN: Community based randomised trial with follow up after 12 months. PARTICIPANTS: 360 patients admitted to three South London Stroke units with stroke or transient ischaemic attack within the past 9 months will be recruited from the wards or outpatients and randomly allocated into two groups. All patients will be visited by the specialist nurse at home at baseline when she will measure their blood pressure and administer a questionnaire. These procedures will be repeated at 12 months follow up by another researcher blind as to whether the patient is in intervention or control group. INTERVENTION: INTERVENTION patients will be given a validated home blood pressure monitor and support from the specialist nurse. Control patients will continue with usual care (blood pressure monitoring by their practice). Main outcome measures in both groups after 12 months: 1. Change in systolic blood pressure.2. Cost effectiveness: Incremental cost of the intervention to the National Health Service and incremental cost per quality adjusted life year gained

Structure Functions and Pair Correlations of the Quark-Gluon Plasma

Recent experiments at RHIC and theoretical considerations indicate that the quark-gluon plasma, present in the fireball of relativistic heavy-ion collisions, might be in a liquid phase. The liquid state can be identified by characteristic correlation and structure functions. Here definitions of the structure functions and pair correlations of the quark-gluon plasma are presented as well as perturbative results. These definitions might be useful for verifying the quark-gluon-plasma liquid in QCD lattice calculations.Comment: 9 pages, 1 figure, revised version (new remark on the coupling parameter on page 2), to be published in Phys. Rev.

Evolutionary optimization of optical antennas

The design of nano-antennas is so far mainly inspired by radio-frequency technology. However, material properties and experimental settings need to be reconsidered at optical frequencies, which entails the need for alternative optimal antenna designs. Here a checkerboard-type, initially random array of gold cubes is subjected to evolutionary optimization. To illustrate the power of the approach we demonstrate that by optimizing the near-field intensity enhancement the evolutionary algorithm finds a new antenna geometry, essentially a split-ring/two-wire antenna hybrid which surpasses by far the performance of a conventional gap antenna by shifting the n=1 split-ring resonance into the optical regime.Comment: Also see Supplementary material, as attached to the main pape

The clinical importance of white matter hyperintensities on brain magnetic resonance imaging: systematic review and meta-analysis

Objectives To review the evidence for an association of white matter hyperintensities with risk of stroke, cognitive decline, dementia, and death

Theory and simulations of rigid polyelectrolytes

We present theoretical and numerical studies on stiff, linear polyelectrolytes within the framework of the cell model. We first review analytical results obtained on a mean-field Poisson-Boltzmann level, and then use molecular dynamics simulations to show, under which circumstances these fail quantitatively and qualitatively. For the hexagonally packed nematic phase of the polyelectrolytes we compute the osmotic coefficient as a function of density. In the presence of multivalent counterions it can become negative, leading to effective attractions. We show that this results from a reduced contribution of the virial part to the pressure. We compute the osmotic coefficient and ionic distribution functions from Poisson-Boltzmann theory with and without a recently proposed correlation correction, and also simulation results for the case of poly(para-phenylene) and compare it to recently obtained experimental data on this stiff polyelectrolyte. We also investigate ion-ion correlations in the strong coupling regime, and compare them to predictions of the recently advocated Wigner crystal theories.Comment: 32 pages, 15 figures, proceedings of the ASTATPHYS-MEX-2001, to be published in Mol. Phy

Influence of molecular temperature on the coherence of fullerenes in a near-field interferometer

We study C70 fullerene matter waves in a Talbot-Lau interferometer as a function of their temperature. While the ideal fringe visibility is observed at moderate molecular temperatures, we find a gradual degradation of the interference contrast if the molecules are heated before entering the interferometer. A method is developed to assess the distribution of the micro-canonical temperatures of the molecules in free flight. This way the heating-dependent reduction of interference contrast can be compared with the predictions of quantum theory. We find that the observed loss of coherence agrees quantitatively with the expected decoherence rate due to the thermal radiation emitted by the hot molecules.Comment: 11 pages, 9 figure

Driven Spin Systems as Quantum Thermodynamic Machines: Fundamental Limits

We show that coupled two level systems like qubits studied in quantum information can be used as a thermodynamic machine. At least three qubits or spins are necessary and arranged in a chain. The system is interfaced between two split baths and the working spin in the middle is externally driven. The machine performs Carnot-type cycles and is able to work as heat pump or engine depending on the temperature difference of the baths $\Delta T$ and the energy differences in the spin system $\Delta E$. It can be shown that the efficiency is a function of $\Delta T$ and $\Delta E$.Comment: 9 pages, 11 figures, accepted for publication in Phys. Rev.

Spherically Symmetric Simulation with Boltzmann Neutrino Transport of Core Collapse and Post-Bounce Evolution of a 15 Solar Mass Star

We present a spherically symmetric, Newtonian core-collapse simulation of a 15 solar mass star with a 1.28 solar mass iron core. The time-, energy-, and angle-dependent transport of electron neutrinos and antineutrinos was treated with a new code which iteratively solves the Boltzmann equation and the equations for neutrino number, energy and momentum to order O(v/c) in the velocity v of the stellar medium. The supernova shock expands to a maximum radius of 350 km instead of only about 240 km as in a comparable calculation with multi-group flux-limited diffusion (MGFLD) by Bruenn, Mezzacappa, & Dineva (1995). This may be explained by stronger neutrino heating due to the more accurate transport in our model. Nevertheless, after 180 ms of expansion the shock finally recedes to a radius around 250 km (compared to about 170 km in the MGFLD run). The effect of an accurate neutrino transport is helpful, but not large enough to cause an explosion of the considered 15 solar mass star. Therefore postshock convection and/or an enhancement of the core neutrino luminosity by convection or reduced neutrino opacities in the neutron star seem necessary for neutrino-driven explosions of such stars. We find an electron fraction Y_e > 0.5 in the neutrino-heated matter, which suggests that the overproduction problem of neutron-rich nuclei with mass numbers around A = 90 in exploding models may be absent when a Boltzmann solver is used for the electron neutrino and antineutrino transport.Comment: 6 pages, LaTex, 3 encapsulated postscript figures, revised and shortened version. Astrophys. J., Letters, accepte